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ABSTRACT 

We use the Spitzer Space Telescope to estimate the dayside thermal emission of the exoplanet 
TrES-3 integrated in the 3.6, 4.5, 5.8, and 8.0 fj,m bandpasses of the Infrared Array Camera (IRAC) 
instrument. We observe two secondary echpses and find relative eclipse depths of 0.00346 ± 0.00035, 
0.00372 ± 0.00054, 0.00449 ± 0.00097, and 0.00475 ± 0.00046, respectively in the 4 IRAC bandpasses. 
We combine our results with the earlier K band measurement of De Mooij et al. (2009), and compare 
them with models of the planetary emission. We find that the planet does not require the presence of 
an inversion layer in the high atmosphere. This is the first very strongly irradiated planet that does 
not have a temperature inversion, which indicates that stellar or planetary characteristics other than 
temperature have an important impact on temperature inversion. De Mooij & Snellen (2009) also 
detected a possible slight offset in the timing of the secondary eclipse in K band. However, based on 
our 4 Spitzer channels, we place a 3a upper limit of \ecos{uj)\ < 0.0056 where e is the planets orbital 
eccentricity and u is the longitude of the periastron. This result strongly indicates that the orbit is 
circular, as expected from tidal circularization theory. 

Subject headings: planetary systems - stars: individual: TrES-3 - techniques: photometric - eclipses - 
infrared 



1. INTRODUCTION 

Among the more than 350 exoplanets known to date, 
transiting hot Jupiters present the first opportunity to 
study and understand the exoplanetary atmospheres. Al- 
though they have masses similar to that of the giant 
planets from the solar system, they orbit extremely close 
to their host star (less than 0.1 AU) and we can ex- 
pect them to be tidally locked to their parent star due 
to a fast tidal synchronisation. Their high equilibrium 
temperature (1000 — 2000 K) and the fact that they 
likely have permanent day/night sides presents interest- 
ing challenges and tests for planetary atmosphere models 
and atmospheric circulation. 

When the secondary eclipse of a transiting system oc- 
curs, it is possible to estimate the flux emitted by the 
day side of the planet relative to the star. The infrared 
Spitzer Space Telescope has been used to detect the flux 
emitted by several exoplanets in its 6 photometric chan- 
nels (3.6, 4.5, 5.8, 8.0, 16.0 and 24.0 /im bandpasses) and 
IRS spectrograph and tabulate the broadband infrared 
spectrum (Charbonneau et al. 2005, 2008; Deming et al. 
2005, 2006, 2007; Grillmair et al. 2007; Richardson et al. 
2007; Harrington et al. 2007; Demory et al. 2007; Knut- 
son et al. 2008; Machalek et al. 2008; Grillmair et al. 
2008). These observations haved pointed to the presence 
of two classes of hot Jupiters: 

One class of planets, including HD 189733b (Deming 
et al. 2006; Grillmair et al. 2007; Charbonneau et al. 
2008; Barman et al. 2008) and TrES-1 (Charbonneau 
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et al. 2005), have emission spectra consistent with stan- 
dard ID cloud-free atmosphere models for these planets 
(Hubeny et al. 2003; Sudarsky et al. 2003; Seager et al. 
2005; Barman et al. 2005; Fortney et al. 2005, 2006a,b; 
Burrows et al. 2005, 2006, 2008). Their infrared spectra 
are dominated by absorption features from CO and H2O. 

The other class, including HD 209458b (Deming ct al. 
2005; Richardson et al. 2007; Burrows et al. 2007a; Knut- 
son et al. 2008), TrES-2 (O'Donovan et al. 2009), XO-lb 
(Machalek et al. 2008), and TrES-4 (Knutson et al. 2009) 
have a temperature inversion between 0.1 — 0.01 bars. 
Water bands that appear in emission instead of absorp- 
tion (Fortney et al. 2006a, 2008; Burrows et al. 2007b, 
2008) are the most likely explanation of their observed 
spectrum. 

The first 6 hot Jupiters planets for which infrared mea- 
surements at two or more wavelengths have been pre- 
sented have shown a connection between their equilib- 
rium temperature and the presence of a temperature in- 
version. It has been proposed (Hubeny et al. 2003; Bur- 
rows et al. 2007b, 2008; Fortney et al. 2008) that these 
two different classes may be linked with TiO and VO 
molecules in the high atmosphere in a gas phase depen- 
dent on its effective temperature, that could lead to a 
temperature inversion from their opacity. However, the 
exoplanet XO-lb (Machalek et al. 2008) does not fit this 
rule, as it is shows evidence of a temperature inversion 
despite levels of irradiation comparable to those of HD 
189733b and TrES-1. 

With a period of only 31 hours (1.30619 days ; Sozzetti 
et al. 2009), TrES-3 (discovered by O'Donovan et al. 
2006) has the shortest period of the known transiting ex- 
oplanets observable with Spitzer, with the exception of 
the recently announced WASP-18b (HeUier et al. 2009). 
The very short period of TrES-3 results in a high level of 
irradiation, with an incident flux of 1.6 10^ergs~^cm~^ . 
Its radius of 1.295 Rjup is larger than predicted by simple 
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models of the structure of these highly irradiated short 
period gas giants. The key to understanding the large 
radius may lie in the composition of the planetary atmo- 
sphere, which dictates the planet's cooling after forma- 
tion and hence its final radius at its current age. Several 
theoretical attempts have been made to propose an addi- 
tional energy source in the planetary interior that would 
combat the planetary contraction after formation. Guil- 
lot et al. (2006) and Burrows et al. (2007b) suggested 
that if the bloated planets have significantly enhanced 
metallicities, the resultant increased planetary opacity 
(and hence reduced contraction rate) could explain the 
large radii. Fressin et al. (2008) confirmed that an evo- 
lution model assuming both a linear correlation between 
the mass of the core (or heavy elements) of giant planets 
and their host star metallacity and an internal energy 
source was likely to reproduce quantitatively the distri- 
bution of masses and radii of the known transiting giant 
planets. 

De Mooij & Snellen (2009) have obtained for TrES-3 
the first measurement of a planetary secondary eclipse 
depth from the ground, in K band using the William 
Herschell telescope (WHT) and the United Kingdom In- 
frared Telescope (UKIRT). They measured the K band 
secondary eclipse depth of 0.241 ± 0.043%. This corre- 
sponds to a day-side brightness temperature at 2.2/xm of 
2040 ± 185 K in _ftr band. They also found the center 
of the secondary eclipse was slightly offset from orbital 
phase 9 = 0.5, of 6*0 = 0.5042±0.0027, indicating that the 
orbit of TrES-3 was perhaps non-circular. If TrES-3 has 
a slightly eccentric orbit, tidal heating from ongoing cir- 
cularization might provide enough energy to explain the 
planet's inflated radius. By measuring the timing of the 
secondary eclipse of TrES-3 in the 4 IRAC bandpasses, 
we will be able to constrain more precisely the planet's 
orbital eccentricity, either confirming or ruling out on- 
going circularization as the explanation for the planet's 
inflated radius. 

2. OBSERVATIONS 

We used the Infrared Array Camera (IRAC ; Fazio et 
al. (2004)) of the Spitzer Space Telescope (Werner et al. 
2004) to observe the secondary eclipse of TrES-3 on UT 
2008 July 18 and July 20, obtaining data at 3.6, 4.5, 5.8, 
and 8.0 fim. We were able to observe it in full array mode 
in all 4 channels for a duration of 5.2 hours We observed 
the target in the IRAC stellar mode, in which the cam- 
era gathers two 10.4 s integrations in the shorter wave- 
lengths channels while gathering a single 30 s integration 
in the longer wavelengths channels. Therefore, we gath- 
ered 1248 images at 3.6, and 4.5 /xmand 624 images at 
5.8 and 8.0 fxm.We describe below our observations in 
two sections, as the InSb detectors used for IRAC chan- 
nels at 3.6 and 4.5 /zm requires a different treatement 
than the Si:As detectors of IRAC 5.8 and 8 /im channels. 
For these observations, we used the "preflash" technique 
(Knutson et al. 2009), in which we pointed the telescope 
towards a bright source before observing TrES-3. This 
was completed in order to reduce the amplitude of the de- 
tector "ramp" at 5.8, and 8.0 /im, effectively pre-loading 
the pixels on which the target star would be pointed. 

2.1. 3.6 and 4-5 jJ-m Observations 
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Fig. 1.— Secondary eclipse of TrES-3 observed on UT 2008 
July 20 at 3.6 and 5.8 ^m, and on UT 2008 Jul. 18 at 4.5 and 
8 ftm. Data are binned in 9.2 minute intervals and normalized to 
one, then offset by a constant for the purposes of this plot. The 
overplotted curves show the best-fit corrections for detector effects 
(see §2.1 and §2.2). 



The contribution of the background to the total flux 
from TrES-3 is low in the 3.6 and 4.5 /zm IRAC band- 
passes, contributing only 0.3% and 0.35%, respectively, 
to the total flux in an aperture with a 5-pixel diame- 
ter centered on the position of the star. We obtain the 
lowest RMS time series using an aperture with a radius 
of 5.0 pixels. We allow the position of our aperture to 
shift with the position of the star in each image. We 
determine the position of the star in each image as the 
position- weighted sum of the flux in a 5-pixel radius disk 
centered on the approximate position of the star. We 
estimate the background in each image from an annulus 
with an inner radius of 12 pixels and an outer radius of 20 
pixels centered on the position of the star. We calculate 
the JD value for each image as the time at mid-expostire 
and apply a correction to convert these JD values to the 
appropriate HJD, taking into account Spitzer's orbital 
position at each point during the observations. 

The most important noise sotirce in the first two IRAC 
bandpasses is due to a well-known intra-pixel sensitiv- 
ity (Reach et al. 2005; Charbonneau et al. 2005, 2008; 
Morales-Calderon et al. 2006; Knutson et al. 2008). 
Fluxes at these two wavelengths show a strong corre- 
lation with the intra-pixel position of the star on the 
detector, at a level comparable to the depth of the sec- 
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Fig. 2.— Secondary eclipse of TrES-3 observed on UT 2008 Jul. 
20 at 3.6 and 5.8 ^m, and on UT 2008 Jul. 18 at 4.5 and 8 ^^m, with 
best-fit eclipse curves overplotted. Data have been normalized to 
remove detector effects (see discussion in §2.1 and §2.2), and binned 
in 9.2 minute intervals, then offset by a constant for the purposes 
of this plot. 

ondary eclipse. We use the following parameters to fit 
the observed fliix as a linear function of the subpixel po- 
sition: 

f = fo * {ci + C2{x - xq) + csiy - yo)) (1) 

where /o is the original flux from the star, / is the 
measured flux, x and y denote the location of the flux- 
weighted ccntroid of the star on the array, xq and yo are 
the coordinates of the center of the pixel containing the 
peak of the star's point spread function, and ci — C3 are 
free parameters in the fit. In the 3.6 /xm channel a;o and 
yo had values of [171.5,175.5], and in the 4.5 /zm chan- 
nel they had values of [167.5,175.5]. We found that the 
position of the star on the array varied by 0.12 pixels in 
X and 0.17 pixels in y during our 3.6 /xm observations. 
During our 4.5 nm observations the position of the star 
varied by 0.11 pixels in x and 0.15 pixels in y. In contrast 
to previous observations of HD 189733 and HD 209458 
in these channels (Kmitson et al. 2008; Charbonneau et 
al. 2008), we find that adding quadratic terms to this 
equation does not improve the fit, likely due to the lower 
SNR of the present observations. In both bandpasses the 
value for the fits is not improved by the addition of 
higher-order terms in x, y and x * y, and the addition 
of these higher-order terms did not significantly alter the 



best-fit values for the best-fit eclipse times and depths. 

After correcting for the intrapixel sensitivity, a trend 
is still visible at 3.6 /zm. Similar observations of TrES-4 
(Knutson et al. 2008) a star with compatible brightness, 
also show the same kind of linear trends at 3.6 /xm and 
4.5 iim, that is likely to be an instrumental effect re- 
lated to the detector or telescope. We correct for this 
effect by fitting the data in both channels with a linear 
function of time. This term is fitted simultaneously with 
the transit curve and the correction for the intrapixel 
sensitivity so that we can accurately characterize the 
additional uncertainty in the depth and timing of the 
eclipse introduced by these corrections. This means that 
at 3.6 /um we are fitting for six parameters, including 
a constant term, a linear function of x position, a lin- 
ear function of y position, a linear function of time, the 
eclipse depth, and the eclipse time. We fit the data using 
a Markov Chain Monte Carlo method (Ford 2005; Winn 
et al. 2007) with 10^ steps, where we set the uncertainty 
on individual points equal to the standard deviation of 
the out-of-transit data after correction for the various 
detector effects. 

Before beginning our fit we do an initial trim within 
our aperture, discarding outliers higher than 3.5 a of the 
local median fiux (defined as the median of a 15 min- 
utes window centered on the data point). We also re- 
move measurements for which the identified position of 
the photocenter x or y deviates more than 3.5 a from 
the same 15-minutes median position. This global trim- 
ing respectivly excludes 6 and 7% of the data points in 
the 3.6 and 4.5 /zm bandpasses. 

Next we carry out the Markov chain fit on the trimmed 
data. We allow both the depth and time of the secondary 
eclipse to vary independently for the eclipses at each of 
the two observed wavelengths, and take the other param- 
eters for the system (planetary and stellar radii, orbital 
period, semi-major axis and inclination) from Sozzetti et 
al. (2009). We calculate our eclipse curve using the equa- 
tions from Mandel & Agol (2002). During each step of 
the chain we exclude outliers greater than either 3.5cr (for 
both the 3.6 and 4.5 /zm fits), as determined using the 
residuals from the model light curve, from our evalua- 
tion of the function. We rescale the value of the 
function to account for the fact that we are varying the 
number of pixels included in the fit. 

After running the chain, we search for the point in the 
chain where the value first falls below the median of 
all the values in the chain (i.e. where the code had 
first found the best-fit solution), and discard all the steps 
up to that point. We take the median of the remaining 
distribution as our best-fit parameter, with errors calcu- 
lated as the symmetric range about the median contain- 
ing 68% of the points in the distribution. The distribu- 
tion of values is very close to symmetric and gaussian for 
the 5 parameters we fitted together (ci , C2 , C3 and tran- 
sit depth and time), and we checked that there were no 
strong correlations between variables. Table 1 states our 
results for the eclipses depths and times and we plot the 
time series in Fig 1. and Fig 2. 

2.2. 5.8 and 8.0 fim Observations 

Previous secondary eclipse studies (e.g. Knutson et al. 
2008) have shown that PSF-fitting can provide a bet- 
ter signal-to-noise ratio at longer wavelengths. At longer 
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wavelengths the flux from the star is smaller and the zo- 
diacal background is larger; we find that the background 
contributes 14% and 16% of the total flux in a 3. -pixel 
aperture at 5.8 and 8.0 /xm, respectively. Because the 
background is higher in these two channels (the median 
background flux is 1.2 MJy/Sr in the 5.8 /xm bandpass 
and 0.6 MJy/Sr in the 8.0 /im bandpass), we used a PSF 
fit to derive the time series in both bandpasses and com- 
pared the results to those from aperture photometry. 

At 5.8 fim and 8.0 /xm respectively, we found that the 
relative scatter in the time series after model fitting from 
the PSF fits was 20% and 25% higher than in the time 
scries from aperture photometry with a 3.0 pixel radius. 
As a result of this increased scatter, which is likely pro- 
duced by discrepancies between the model PSF and the 
observed PSF, we conclude that aperture photometry is 
also preferable in these 2 channels. We compare the time 
series using apertures ranging from 3 — 4.5 pixels and find 
consistent results in all cases, but with a scatter that in- 
creases with the radius of the photometric aperture. 

The 5.8 and 8 /xm bandpasses of IRAC camera are 
known to be affected by a "detector ramp" (Kmitson et 
al. 2008, 2009; Charbonneau et al. 2008) that causes the 
effective gain (and thus the measured flux) in individ- 
ual pixels to increase. The size of this effect depends on 
the illumination level of the individual pixel. Pixels with 
high illumination (> 250M Jy / Sr in the 8 /im channel) 
will converge to a constant value within the first hour of 
observations, whereas lower-illumination pixels will show 
a linear increase in the measured flux over time with a 
slope that varies inversely with the logarithm of the illu- 
mination level. Preflashing the detector array by point- 
ing a bright source prior to conducting the science obser- 
vations can reduce the amplitude of this ramp by doing 
a pre-load of the pixels on which the target star will 
be pointed. In our observation of TrES-3, we targeted 
a bright star as a preflash source and this yielded one 
of the smallest ramp effect ever recorded using 5.8 and 
8 /im observations. Previous observations (e.g. Knutson 
et al. (2009)) have shown that the ramp is well described 
as following an asymptotic shape, with a steeper rise in 
the first 30 minutes of observations. We correct for this 
effect by fitting our time series in both bandpasses with 
the following function: 

/ = /o * (ci + c^lnidt)) (2) 

where /o is the original flux from the star, / is the 
mcasiired flux, and dt is the elapsed time in days since 
the start of the observations. We fit both Eq. 2 and the 
transit curve to the data simultaneously using a Markov 
Chain Monte Carlo method as described in §2.1. As be- 
fore, the distribution of values was very close to sym- 
metric in all cases, and there were no strong correlations 
between the variables. Best-fit eclipse depths and times 
from these fits are given in Table 1, and the time se- 
ries before and after correcting for detector effects are 
shown in Figures 1 and 2, respectively. As a check we 
repeated these fits adding a quadratic term of ln{dt) in 
Eq. 2, and found that the value of the function for our 
best-fit solution were similar at both 5.8 and 8.0 /xm to 
their previous values. In order to estimate error bars for 
the measured depth and time, we take the median of the 
remaining distribution as our best-fit parameter, with er- 
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Fig. 3. — Probability distribution for the eclipse depth from 
a Markov Chain Monte Carlo fit to the IRAC 8 /im data. The 
long-dashed, short-dashed and dot-dashed a lines indicate the la, 
2(T and 3a limits on the eclipse depth, respectively, which were 
calculated by integrating over this distribution. 

rors calculated as the symmetric range about the median 
containing 68% of the points in the distribution. Figure 
3 shows the representative example of the histogram of 
the probability distribution for the eclipse depth at 8.0 
/xm. 

3. DISCUSSION 

We determine the best-fit eclipse times for the four 
secondary eclipses observed using IRAC by taking the 
weighted average of the best-fit eclipse times in each 
bandpass. Using this method, we find that the eclipse is 
shifted by LO^i J minutes later than the expected time 
based on the ephemeris from Sozzetti et al. (2009). 

Our estimate for the best-fit timing offset translates to 
a constraint on the orbital eccentricity e and the argu- 
ment of pericentcr uj of ecos(a;) = 0.00084^Q QQQg; the 3cr 
upper limit on its absolute value is \ecos{u)\ < 0.0056. 
We selected these limits because we are interested in con- 
straining the magnitude of e rather than the sign of the 
cos{lu) term. This upper limit means that unless the 
longitude of periastron lj is close to 90° or 270°, we can 
rule out tidal heating from ongoing orbital circulariza- 
tion (Bodenheimer et al. 2001; Liu et al. 2008) as an 
explanation for the infiated radius of TrES-3. Winn et 
al. (2008) have placed a constraint on the albedo as- 
suming a circularized orbit. Thus, our results serve to 
validate this assumption, and we proceed to adopt their 
stated upper limits in our analysis. They placed a 99%- 
confidence upper limits on the planct-to-star fiux ratio 
of 2.4 X 10"'', 5.0 X lO"'', and 8.6 x lO"'' in the i, z, and 
R bands respectively. 

We then compare the secondary eclipse depths in the 
four IRAC bandpasses and the previously reported K 
band value to the predictions from atmosphere models 
for this planet (see Fig. 4 and 5). We employ the formal- 
ism described in Burrows et al. (2007a, 2008). Wc use a 
stellar atmosphere model (Kurucz 1979, 1994, 2005) with 
an effective temperature of 5650 K and a planet-star ra- 
dius ratio of 0.1654 based on the measures of Sozzetti 
et al. (2009). We calculate the emergent spectrum at 
secondary eclipse for a pair of free parameters, P„ and 
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TABLE 1 
Best-Fit Eclipse Depths and Times 



A (/im) 


Eclipse Depth 


Center of Transit (HJD) 


O-C (min.)'' 


3.6 


0.356 ± 0.035% 


2454668.5447 ±0.0020 


-2.0 ± 2.9 


4.5 


0.372 ± 0.054% 


2454665.9343 ± 0.0027 


0.9 ±3.9 


5.8 


0.449 ± 0.097% 


2454668.5498 ± 0.0042 


5.4 ±6.0 


8.0 


0.475 ± 0.046% 


2454665.9365 ± 0.0021 


4.0 ±3.0 



^Observed minus calculated transit times, where the expected transit times 
are calculated using the ephemeris from Mandushev et al. (2007) and assum- 
ing zero eccentricity. 
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Fig. 4. — Day-side planet-star flux ratios for TrES-3 as determined from measurements of the secondary eclipse depth in the four IRAC 
bandpasses (blue circles) and in the K band UKIRT detection. We also show the three 99 %-confidence upper limits obtained by Winn 
et al. 2008 in the R, i, and z bands. The black line corresponds to a default model with no temperature inversion with a redistribution 
parameter P„ = 0.5, which describes the case where the incident energy is fully redistributed across the entire surface of the planet. The 
red and green lines correspond to models with an additional optical absorber at high altitudes (parameterized as Kextra), which produces a 
thermal inversion around pressures of 0.001 bar (Burrows et al. 2007a, 2008). Squares show the values for these models after integrating over 
the Spitzer bandpasses. The eclipse depths in the 3.6, 4.5 and 5.8 /im bandpasses can be matched by both models with relatively efficient 
day-night circulation P„ = 0.3 and different levels of additional opacity (red model Kextra = 0.1 cm^/g and green model Kextra = 0.01). 
However, the planet-star ffux ratio at 8.0 /im argues strongly against the presence of a temperature inversion, as the green model with a 
low Kextra = 0.01 cm^/g provides the best match at this wavelength ; the K-band UKIRT detection also slightly favors this scenario. 



Kextra- Pn IS the dimensioiiless redistribution parame- 
ter that accounts for the coohng of the dayside and the 
warming of the nightside by zonal winds near an opti- 
cal depth of order unity, ranging form to 0.5. It is a 
measure of the efficiency of heat redistribution by super- 
rotational hydrodynamic flows. As the value of P„ is 
increased, the day side becomes cooler and the emer- 
gent planetary flux at superior conjunction becomes cor- 
respondingly small. Kextra is thc absorptive opacity in 
the optical at altitude (here in cm^/g) and our best fit 
requires only a very low opacity that does not create a 



temperature inversion. 

Fig. 4 and 5 show three models with different val- 
ues for Pn and Kextra- Thc standard non-inverted model 
{Kextra = cm^/g) is clcarly inconsistent with the ob- 
served fluxes from TrES-3 at wavelengths longer than 
3 /im. It is possible to match the observed 3.6 /xm flux 
with this model by reducing the relative fraction of the 
incident energy that is redistributed to the planet's night 
side, thus increasing the day-side temperature and cor- 
responding fluxes. We plotted the best-fit model, that 
involves no thermal inversion {Kextra — 0.01 cm^/g), and 
an example of an inverted {Kextra — 0.1 cm'^/g) model. 
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TABLE 2 

VALUES FOR DIFFERENT 
ATMOSPHERIC MODELS 



f^extra 


Pn 






0.00 


0.1 


1.317 


1.105 


0.01 


0.1 


3.895 


1.289 


0.1 


0.1 


41.10 


13.70 


0.00 


0.5 


40.77 


13.59 


0.01 


0.5 


24.60 


8.210 


0.1 


0.5 


38.43 


12.76 


0.00 


0.3 


6.346 


2.115 


0.01 


0.3 


1.734 


0.578 


0.025 


0.3 


2.645 


0.881 


0.05 


0.3 


9.400 


3.132 


0.1 


0.3 


25.53 


8.511 



^X^/i^ is the reduced with 
n — 5 parameters —2 degrees of 
freedom. 
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Fig. 5. — Day-side pressure-temperature profiles for the three 
models plotted in Fig. 4. The temperature increases at low pres- 
sures as a function of the amount of the opacity of the Kextra 
absorber. The fraction of energy redistributed to the night side of 
the planet Pn decreases the temperature at lower atmospheric lev- 
els (0.01 - 0.1 bars). A drop in temperature (described in Burrows 
et al. 2007b, 2008) occurs as the day-night circulation is turned 
up. 

Although, for each of the four individual planet-star flux 
ratios we are able to find a model with a temperature 
inversion that fits well, we are unable to find a single 
model that fits all four data points simultaneously. In 
the (Kextra = 0.1 cm^/g) sceuario involving a tempera- 
ture inversion, the 8.0 /xm flux can not be reproduced 
well, and we concluded that the model with no inversion 
provides an overall better fit. The K band measurement 
also strengthens the conclusion that the model with no 
inversion provides a better global fit to the data. We 
compared the values for 5 data points and 2 free pa- 
rameters for several models including the three plotted 
models: table 2 shows the values we obtain these 
different models. Only non-inverted models provide a 
good fit to the 5 data points, with the plotted model 
[i^extra — O.OI cm^/g and P„ = 0.3) showing the best 
solution. 

4. CONCLUSIONS 



We have detected the TrES-3 secondary eclipse in the 4 
bandpasses of the IRAC instrument. These observations 
at 3.6, 4.5, 5.8, and 8.0 /xm, combined with the K band 
measuement from De Mooij & Snellen (2009), reveal that 
this planet does not show a thermal inversion similar to 
the one observed for HD 209458b (Knutson et al. 2008; 
Burrows et al. 2007a), TrES-2 (O'Donovan et al. 2009) 
and TrES-4 (Knutson et al. 2009). The best overall fit 
involves an efficient day-night circulation (P„ = 0.3) and 
a very low additional opacity {Kextra = 0.01 cm^/g). The 
scenarios presented by Fortney et al. (2008) would pre- 
dict that gas phase TiO or VO at high altitude would re- 
sult in a temperature inversion for this highly-irradiated 
temperature planet, as it is warmer than IID209458, 
TrES-2 and TrES-4, which are inverted. The fact that 
our results strongly favor a scenario without any temper- 
ature inversion shows that the distinction is not simply 
due to the level of irradiation for the separation between 
these two kinds of close-in giant planets. Spiegel et al. 
(2009) have shown that the TiO-VO hypothesis was un- 
likely unless there is significant mixing in the atmosphere, 
and TrES-3 could be among the planets were the mixing 
is insufficient. 

With an upper limit of |ecos(a;)| < 0.0056 for the or- 
bital eccentricity, we can rule out tidal heating from on- 
going orbital circularization at the level required by Liu 
et al. (2008) in order to explain TrES-3's inflated ra- 
dius. This is also the flrst exoplanet in this range of 
high irradiation level not to show a temperature inver- 
sion. This interesting and unexpected case emphasizes 
the importance of gathering more hot Jupiter infrared 
emission measurements, in order to study the correla- 
tions between temperature inversion and system char- 
acteristics. A large sample of multi- wavelength infrared 
measurements from many different exoplanets will be re- 
quired to understand the origin of these temperature in- 
versions. Fortunately, this will be accessible during the 
Warm phase of the Spitzer mission, as both the 3.6 and 
4.5 fim channels will continue to function at full sensitiv- 
ity. Thirty-four known transiting exoplanets known to 
date are bright enough for the two Spitzer observations 
at 3.6 and 4.5 /im to assess or rule out if a temperature 
inversion occurs in their high atmosphere. Observations 
of the secondary eclipse in these two channels should be 
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sufficient to distinguish between planets with and with- 
out temperature inversions in most cases, ahhough the 
case of TrES-3 shows that the additional information at 
longer wavelengths can be the discriminating factor be- 
tween these two possible scenarios. The ground-based 
high precision detection of thermal emission is also help- 
ful. The global set of observations, combining detections 
in the 4 IRAC bandpasses obtained during the cryogenic 
Spitzer mission for 15 planets, and the observations of 19 
other planets during the warm phase of the mission by 
the Spitzer Exploratory Science Program 60021 (Princi- 
pal Investogator H. Knutson) will allow the study of cor- 
relations of the temperature inversion phenomenon with 



stellar metallicities, planet size and mass, levels of irra- 
diation, surface gravities, and orbital periods. Studying 
the connection between the presence of a temperature 
inversion and these characteristics will give us a better 
understanding of planetary atmospheres under strong ir- 
radiation. 



This work is based on observations made with the 
Spitzer Space Telescope, which is operated by the Jet 
Propulsion Laboratory, California Institute of Technol- 
ogy, under contract to NASA. 
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